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Abstract 

The non-perturbative Wegner-Houghton renormahzation group is analyzed by the local po- 
tential approximation in O(A^) scalar theories in d-dimensions (3 < c? < 4). The leading 
critical exponents v are calculated in order to investigate the effectiveness of the local poten- 
tial approximation by comparing them with the other non-perturbative methods. We show 
analytically that the local potential approximation gives the exact exponents up to 0(e) 
in e-expansion and the leading in 1/A^-expansion. We claim that this approximation offers 
fairly accurate results in the whole range of the parameter space of N and d. It is a great 
advantage of our method that no diverging expansions appear in the procedure. 



1 Introduction 

The non-perturbative phenomena of the quantum field theories have been fascinating 
many physicists. There are only a limited number of tools to attack such problems, for 
example, the Monte Carlo simulations in the lattice field theories, the Schwinger-Dyson 
equations, e-expansion, 1/A^-expansion etc. In this article we focus our attention on the 
Wilson renormalization group (RG) among these approaches |l|]. The Wilson RG equations 
are given by the form of the functional differential equations for the so-called Wilsonian 
effective action defined in the Euclidean space. Therefore it is inevitable to approximate 
them for the practical calculations. We usually expand the effective action in terms of the 
number of derivatives included in the general operators and solve the Wilson RG equation 
within the subspace up to some finite number of derivatives. In the first order of this 
approximation, any derivative couplings are dropped except for the kinetic terms. This is 
called the local potential approximation (LPA). Then the functional differential equation 
for the effective action is reduced to a non-linear partial differential equation for the local 
potential. Instead of analyzing this partial differential equation directly we may also expand 
the local potential with respect to the fields and truncate the series at some finite order. 
By this approximation we may solve the coupled differential equations for the expansion 
coefficients, that is, the coupling constants. These equations are much easier to solve than 
the original partial differential equation. Actually this series expansion is found to converge 
remarkably fast. 

It is important to see whether the LPA offers us sufficiently good results before performing 
the higher order calculations in the derivative expansion. For this purpose we take the 
Wegner-Houghton (W-H) equation 0] among several formulations]^, ^ and compare the LPA 
of this with e-expansion as well as 1/A^-expansion in O(A^) scalar field theories. 

The contents of this article are the following. In section 2 we briefiy review the W- 
H equation for O(A^) scalar field theories and the LPA. We series expand the potential 
and calculate the correlation length critical exponent. We compare our results with the 
e-expansion in section 3 and the 1/A^-expansion in section 4. The accuracy of the LPA 
is discussed in comparison with the Borel resummation of the e-expansion and other non- 
perturbative calculations. We also argue the property of the LPA in relation to the 
expansion. The section 5 is devoted to the summary and discussions. 



2 The Wegner-Houghton equation and the local po- 
tential approximation 

First we derive the W-H equations for the 0{N) scalar theories^]. The starting point is 
the Euclidean path integral defined with momentum cutoff A{t) = e^*A; 



where 5'eff is the Wilsonian effective action. The RG differential equations give the response 
of Ses under the infinitesimal change of the cutoff A with keeping the partition function 
Z unaltered. The difference in S^s induced by the change of the cutoff is determined by 
integrating the "shell mode" with the momenta between A{t) and A{t + 6t). It should be 
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noted that this integration is reduced to the Gaussian one for infinitesimally small St, and is 
exactly carried out. Besides we rescale the momentum and the fields by A, since the change 
of the dimensionless quantities are of our interest. Thus we deduce the W-H equation for 
O(A^) scalar field theories as 

dS^s 1 f d'^p / S'^S.s \ 

In- 



dt 26tJe~st^\p\<i{27rY \ 6<P{p)6<j){-p) J 

_ J_ f f d'^p d^q [ 6S,s ( S^S,s y SS,s \ 

25t J ye-«<b|,|,|<i (27r)^ (27r)^ \ d^ip) J 

Here 0" denote iV-component scalar fields and 7]/ 2 is the anomalous dimension of these 
fields. 

The effective action S'efr is given in the derivative expansion by 

^eff = / d'x {V{p) + ^Z,ip){d,rf + \z2{pWd'r + } , (3) 

where we introduced the O(A^) invariant operator p = If we substitute this into 

Eq.(||), the W-H equation turns out to be the coupled partial differential equations for V{p), 
Zi{p) (z = 1, 2, ...). In the LPA we restrict the S'efr to 



'eff 



Jd'^x[v{p) + l{d,r)'Y (4) 



Note that the kinetic term is not renormalized in the LPA, which means that the anomalous 
dimension automatically vanishes. Then the W-H equation in the LPA is written down as 

^ = ^ [ln(l + V' + 2pV") + {N- l)ln(l + V')] + d-V + (2 - d)pV\ (5) 

where the prime denotes the derivative with respect to p and is the d dimensional angular 
integral, = n-'^/^2^-'^/T{d/2). 

Next we are going to examine this partial differential equation by expanding the potential 
into a power series and by truncating it at some order. Such analysis based on the Taylor se- 
ries will be useful when comparing the LPA with the e-expansion and also the 1/A^-expansion 
as we will see later. If the truncation method gives the same results irrespective to the order 
of the truncation, we consider that they are the LPA results. 

The naive expansion may be to expand around the origin (Fixed Scheme), 

Vip) = E ^P-. (6) 

Also it will be natural to expand around the minimum of the potential bi(t) (Comoving 
Scheme) P] , 

v{p) = bo{t) + j:^-^{p-h{t)r. (7) 

m=2 
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These two schemes give the RG flows projected onto the different M-dimensional subspaces. 
In Fig.l, we show the leading critical exponents v (the correlation length exponents) esti- 
mated by these two truncation schemes for A^=l theory in (i=3. It is seen that the Comoving 
Scheme converges very rapidly to a value 0.68956(1), which should be compared with 0.687(1) 
obtained by the analysis of the partial differential equation (5)0. Contrary to this good 
convergence, the exponents estimated by the Fixed Scheme seem to remain oscillating. The 
reason of this oscillating behavior has been discussed by Morris^]. Anyhow this convergence 
indicates that the relevant operator can be described sufficiently well within quite a few 
dimensional subspace given by Eq.(7)[H]. 
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Figure 1: The critical exponents for the A^=l scalar theory in d—Z calculated by the Fixed Scheme and the 
Comoving Scheme as a function of the order M of the truncation. 



Now let US consider general coordinates introduced to describe the RG flows. Suppose 
M. is the theory space of infinite dimension and {a*} is a generic coordinate system (a set 
of coupling constants) of M.. The operator djdt defines a vector field in M. and is given in 
the coordinate system {a*} by 

where we define the generalized beta functions /5*(a) = da^/dt. This means that the beta 
functions are the components of the vector field d/dt. For any function S{a^) defined in M., 
dS/dt gives a tangent vector in the space of {S}; 



dS_ 
It 



(9) 



where we also introduced the base vectors in this space, = dS/da\ In practice the RG 
equation gives us dS/dt, where S stands for the Wilsonian effective action, and we may 
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obtain the beta functions by expanding dS/dt by the base vectors ^j. Here it should be 
noted that the base vectors depend on a*, that is, the position in generically. We may 
call such coordinates system a "comoving frame" . The expansion at the potential minimum 
given in Eq.(|^) is a typical example of this comoving frame, since the base vectors given 

explicitly by ^1 = — 62 — &3(p — &i) , = p — ^i, ■Cs = (P"~^i)^5 depend on the coordinates 

bi indeed. 

Suppose we may choose the coordinates {a*} so that the matrix fi* = df3^/daj forms a 
lower triangle, 

^(a) = |^(a) = for t < j, (10) 

then the components a*(t) of the flow (curve) can be exactly evaluated within the finite 
dimensional subspaces of A4 spanned by (a^, a^, ■ ■ ■ a*). Here let us call such a special 
coordinate system the "perfect coordinates". Generally it would be difficult to find such 
coordinates as much as to solve the RG equation exactly. However in the large limit the 
expansion around the potential minimum will be found to give us an example of the "perfect 
coordinates" . Actually we can solve the RG equation in large N limit exactly in the every 
order of the truncation. 



3 The comparison with the e-expansion. 

The critical exponents have been calculated in powers of e = A — d hj Wilson and 
Fisher |l|. The correlation length exponent u is found to be 



N + 2 ^_^{N + 2){N'^ + 23N + 60) ^2 
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2 ^ 4(iV + 8) 



0(e' 



]{N + 8f ^ — vw (11) 

In this section we consider to compare the exponent obtained by the W-H equation in the 
LPA with this result. 

We examine Eq.(5) in d = 4 — e. In the Fixed Scheme defined by Eq.(6) the first three 
beta functions are given by 
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From these equations the non-trivial fixed point solution is found to be 



+ 2(e-l)o3. 
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(13) 
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d 


4.0 


3.9 


3.8 


3.7 


3.6 


3.5 


3.4 


3.3 


3.2 


3.1 


3.0 


V 


0.500 


0.509 


0.519 


0.531 


0.544 


0.560 


0.577 


0.598 


0.622 


0.652 


0.6896 



Table 1: The critical exponents for the A^=l scalar theory in d-dimension. 




Figure 2: The critical exponents for the A''=l scalar theory as a function of e. The estimations by the 
W-H equation in the LPA and by e-expansion up to 0(e') are shown. The bold solid line is by the Borel 
resummation of the e-expansion. 



in powers of e, where A/^ = l/87r^, A'^ = (1 — 7b + ln47r)/167r^ and •je = 0.5772.... denotes 
the Euler constant. The exponent u is given by the positive eigenvalue of the matrix f2* = 
d(3i/daj at the non-trivial fixed point. Thus the exponent estimated in the LPA turns out 
to be 

1 N + 2 {N + 2){N^ + 38N + 96) 2 3^ 

Comparing with Eq.(ll) the LPA exponent is exact in O(e)0. The deviation starts in O(e^), 
which is due to the ignorance of the renormalization of the derivative couplings. 

We list the exponents calculated by the W-H equation in Table 1, and show them in 



Fig. 2 as well as the exponents obtained by e-expansion for A^=l[|l^, |Tl[. The e-expansion 
itself produces an asymptotic series, which means the series converges up to a certain order 
but eventually turns to diverge. Therefore the e-expansion would not be meaningful, unless 
this asymptotic series is summed up to a finite value. Actually not only the e-expansion but 
also the ordinary perturbation with respect to the coupling constants has been known to 
generate an asymptotic series. Since the late 70's the Borel summation technique has been 
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found to be quite powerful to sum up such asymptotic series and has succeeded in offering 



the rehable values for the critical exponents ||T2|, |10| |13| . The critical exponent v calculated 



by the Borel summation technique a ia [|1^, |T3| as a function of e is also drawn in Fig.2. This 
may well be regarded as the accurate result for the exponent with which we can compare our 
RG calculations. Then our method gives comparatively good values for the exponent, if we 
take into consideration that the LPA is just the first order approximation in the derivative 
expansion. 

Of course it will be extremely important to investigate whether the Wilson RG calculation 
in the higher order approximation than the LPA really gives the results converging order by 
order. However it should be noted that the Wilson RG seems to give much easier methods 
to handle than the e-expansion and also than the perturbative expansion at least in order to 
obtain moderately reliable results. Because the higher order of such expansions necessitates 
quite complicated calculations. Moreover it is essential to know the large order behavior 
of the asymptotic series in performing the resummation. Besides it would be necessary 
to tune some free parameters introduced to define the Borel transformation so as to reach 
the accuracy shown in Ref . pUj, |TB[. Thus the procedure of the Borel summation is rather 
complicated. On the other hand no such expansions leading to asymptotic series are carried 
out in our calculation of the Wilson RG equation. Therefore we may say that our RG method 
offers directly the results which are supposed to be sums of asymptotic series appearing in 
other methods. Thus this is one of the characteristic and advantageous features of the Wilson 
RG approach. 



4 The comparison with the 1/A^-expansion. 



It has been recognized p|, |T4[ that the LPA for the A^- vector model gives the exact effective 



potential in large A^ limit. To see this, we expand the effective action as S'efr = Kfr + 
S'gg + S'gg + where the index denotes the number of the operator 0^ with no n- vanishing 
momentum. The beta-functional of V^s depends only on {Kfr} in large A^ limit 0. That is, 
Kff can be seen as a sort of the first component of the perfect coordinates introduced in 
Eq.(lO). Therefore we may solve the RG equation exactly in the subspace irrespectively 
to other subspaces. The RG equation for this subspace turns out to be the W-H equation 
d^). This feature is also expected from the fact that the renormalization of the derivative 
couplings is required only in 0(1/A^)[|r^, which is discarded in the LPA. 



Here we also show this remarkable property by deducing the RG equation directly from 
the large A^ effective potential In order to define the large A^ limit RG equation, we 
introduce x = p/N and F{x,t) = dV{p)/dp, and rewrite Eq.(5) as 

dF , , Ad { 3F' + xF" , , F' \ . ,^ , 

— ix^t) = —\ + A^ - 1 ^ + 2F + 2 - (£)xF'. 15 

dt^ ' ' 2A^ \ 1 + 2F + xF' ^ ^ 1 + F J ^ ' ^ ' 

In large A^ limit this equation is reduced to 

dW A P' 

— (x, t) = + 2F + (2 - d)xF'. (16) 

Now we are going to show this equation can be derived directly from the gap-equation 
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in large limit as well. First introduce two auxiliary fields x ^i-nd p, 

Z = j VrVxVpexp |- J d^x \^{dr? + x[l ir? - iVp) + NV{p) 



(17) 



so that we may treat the generic form of the potential. Then the effective potential V^s, to 
say strictly, the constraint effective potential |T^, for A^- vector model may be written as 

exp (-Weff (p)) = J VxVpexp!^-N J d^x [x (p - p) + V{p)] - ytrln(-92 + x)} , (18) 

where p denotes (j)^/2N. In large N limit, the path integral of these auxiliary fields is 
evaluated by the saddle point method. Then the effective potential V^s will be given by 
solving the following coupled equations; 

Kfr[p] = x{p-p) + V{p) + \j^^^Hk' + x\ (19) 
I [ d'^k I 

X - (21) 
dp 

Here the cut-off A in the momentum integration is introduced to derive the RG equation for 
the effective potential Ves. By considering the infinitesimal change of the cut-off, the RG 
equation is found to be 



dA dA 



-^A'''\n{A' + x), (22) 



where the saddle point equations dVcs /dp = and dVes /dx = are used. Noting that x is 
given in terms of Vcs as x = dVcs/dp, we eventually obtain the RG equation for the effective 
potential Ves as 

^[p;A] = -Y^^-Mn(A^ + 0, (23) 

which is found to be identical to the W-H equation in the LPA (p!6D, after rescaling of the 
variables. 

Next we solve this RG equation by projecting it on much smaller subspaces. Actually 
if we expand the effective potential in powers of p — bi{t), where bi is the minimum of the 
potential, as is given by Eq.(|^), then = dbk/dt {k = 1, 2, 3, ...) are given by 

A = -^ + (t/-2)6i, 

P2 = -^bl + {4-d)b2, (24) 
P3 = -Y^3b2b3 + 2bl) + (6 - 2d)b3, etc. 



Here note that, contrary to Eq. (|T^) , /^^ depends only on the couplings bi of i < k. Namely 
this "comoving frame" is found to be the "perfect coordinates" in large A^ limit. This means 
that the flow of bi can be exactly determined independently of the higher order couplings bk 
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oi k > i. In other words the projection to each order of truncated subspace in the comoving 
frame is always exact. Therefore the coordinates of the fixed point are also exactly solved in 
each order of the truncation and are found to be 



Ad 


2{d- 


2)' 


2(4- 


d) 


Ad 




4(4- 


df 



(25) 



etc. 



A%d-6y 

In this perfect coordinates the matrix = djSi/dbj is characterized by a lower triangular 
form. At the non-trivial fixed point, Q is given by 
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n 



d-2 








d-A 






d-6 
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d-8 










v 



(26) 



where * denotes some functions of d. The eigenvalues of this triangular matrix are nothing 
but the diagonal elements, d — 2m, (m = 1, 2, 3, ..). Therefore the leading critical exponent 
in large N limit is exactly derived as 

^ (27) 



d-2 



Now we study the 1/A^ dependence. In d=3, the critical exponent u has been calculated 
up to 0{l/N^Wi 
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In Table 2, we list the critical exponents calculated by our RG method, and show them in 
Fig. 3 as well as those obtained by 1/A^-expansion and by other methods. The RG equation 
in the LPA (5) is fairly effective in the whole range of N, while 1/A^-expansion results show 
up its divergent nature and there has been no way to sum up the series. 



N 


1 


2 


3 


4 


5 


10 


20 


50 


100 


V 


0.6896 


0.767 


0.826 


0.865 


0.891 


0.946 


0.974 


0.990 


0.995 



Tabic 2. The critical exponents of the 0(iV)-synimetric scalar theories in (i=3. 

Indeed the LPA seems to be inferior to 1/A^-expansion for a large but finite A^. This is 
because the estimation of v by the RG deviates considerably from the exact value in 0(1/A^). 
Of course it is well expected that the next-to-leading calculation in the derivative expansion 
will become fairly accurate in this region. However such study has not been carried out 
yet. Here we numerically evaluate the coefficient of 0(1/A^) for the exponent v calculated 
by the W-H equation for 3 < < 4. The results are shown in Fig. 4 with the coefficient of 
1/A^ in Eq.(27). The difference is getting smaller and smaller as the dimension d approaches 
to 4, because the effects of the derivative couplings decrease monotonically and eventually 
vanishes in ci = 4 where the LPA is exact for any A^. 
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Figure 3: The critical exponent of the 0(Af)-symnietric scalar theories in d = 3 as a function of The 
results by the W-H equation in the LPA are shown with the results of 0(l/iV) and 0(l/iV^) in l/A^-expansion. 
The square points show the present best estimates summarized in Ref.[n8[. 
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Figure 4: The dimensional dependence of the coefficient of in the W-H equation and in 1/A^-expansion. 
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5 Summary and discussions 



In this article it has been shown explicitly that the LPA of the W-H equation becomes 
exact in 0(e) and in large limit. This is because the higher order corrections involving 
the derivative interactions are not generated in these order calculations. The deviations 
appearing in the next orders in these expansions are also studied. In Fig. 5 the global behavior 
of the exponents obtained by the W-H equation in the LPA as well as the accurate ones are 
shown so as to summarize our numerical results. It should be noted here that the RG method 
is rather effective even in the LPA over the whole parameter region of A^ and d. Thus we 
can comprehend the gross properties of the exponent by means of our simple Wilson RG 
method. 




• LPA W-H 

X World best estimates 
^ Exact 

— 0(80, 0(1/N2) 



Figure 5: The variation of the exponent with respect to N and d. The fifth order of e-expansion and the 
second order of 1/A^-expansion are plotted for comparison. The world best estimates are taken from Ref.[Q 



The perturbative expansion as well as e-expansion leads to the asymptotic series. There- 
fore such expansions cannot be applied directly at the region where the expansion parameters 
are not small enough. It is necessary to sum up these series through the Borel transformation 
in order to obtain quantitatively reliable results. Moreover we have to know the large order 
behavior of the asymptotic series to carry out the summation. In contrast to these expansion 
schemes, it is remarkable that we can estimate the exponents by solving the W-H equation 
without resorting to the complicated manipulations of the resummation. 

Needless to say it is quite important to see whether the higher order calculations in the 
derivative expansion generate converging results or not. So far the Wilson RG equations 
have been investigated up to the next-to-leading order of the derivative expansion ||19||. By 



these analyses the exponents have been found to become closer to the expected values in 
0((9^) than in the LPA. It should be noted that the sharp cutoff scheme applied to derive the 
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W-H equation leads to the difficulties of non-analyticity when the derivative couplings are 
incorporated [^]. Therefore the formulations defined by the smooth cutoff Q have been 
used in the next order calculations]^. 

The derivative expansion is not an "expansion" . It is just enlarging the functional sub- 
space step by step in which the Wilson RG equation is solved. Therefore it is much more 
like increasing the total lattice size in the Monte Carlo simulations. Nobody may expect 
any divergent behavior of the physical quantities when increasing the lattice size. Rather we 
should expect convergence, or oscillation at worst. We have already met an example of this 
type of convergence in Fig.l, where we enlarge the dimension of the functional subspace one 
by one and we get the strongly converging results (or oscillation, depending on the choice of 
the truncation scheme). Thus our method of the Wilson RG equation may not suffer from 
any divergent series. 

Lastly let us discuss one of the applications to particle physics. The perturbative beta 
functions have been known to show the behavior of the asymptotic series as well. Therefore 
we cannot assert even the existence of the non-trivial fixed point by using the perturbative 
beta functions. Indeed the Borel summation technique has been succeeded in two and three 
dimensions. In four dimensions, however, even the Borel summability for the scalar theories 
has not been clarified yet. Therefore the other non-perturbative methods are desirable in 
four dimensions. For instance the argument of the triviality mass bound for Higgs boson 
and for top quark relies entirely on the non-existence of the non-trivial UV fixed point in 
the standard model. Actually several works on the triviality bound by means of the Wilson 



RG equations have been already done|^. However more minute and realistic investigations 
would be desired [p2| . 

K-I. A. is supported in part by the Grand-in-Aid for Scientific Research (^^04640297) 
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